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^ I Abstract 

' In this paper we study a 1+1 dimensional relativistic parton model for the 

structure of baryons. The quarks and anti-quarks interact through a linear po- 

' tential. We obtain an analytic formula for the isospin averaged valence quark 

^ ■ ..... .1 

. distribution in the chiral and large Nc limits. The leading — and non-zero 

' current quark mass corrections are estimated. Then we extend this model to 

^ i' include 'sea' and anti-quarks. We find that the anti-quark content is small at 

P> i' a low value of Q^. Using these distributions as initial conditions for evo- 

(~| , lution, we compare with experimental measurements of the structure function 

^ ■ xF^^x, Q^) and find reasonable agreement. The only parameters we can adjust 

' are the fraction of baryon momentum carried by valence quarks and the initial 

H ' 

\ scale Qq. 
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1 Introduction 



We present a variational parton model description for the structure of baryons 
as measured in Deep Inelastic Scattering. This model enables us to calculate 
the xb dependence of the structure function xF^ at an initial value of Qq. We 
compare this prediction with experimental measurements by CCFR and CDHS 
collaborations. 

In Deep Inelastic Scattering, the longitudinal momenta of the partons dom- 
inates their transverse momenta. Indeed, as pointed out by Altarelli, Parisi and 
others ^ , a perturbative treatment of transverse momenta, with an upper cut- 
off Q, leads to the same scaling violations as predicted by the Operator Product 
Expansion in the leading logarithmic approximation. By the uncertainty prin- 
ciple, the virtual photon momentum Q, is a measure of the size of transverse 
momenta being probed. 

While the dependence (for sufhciently large Q^) of the structure functions 
is well understood IT]|2|, the xb dependence is harder to understand since it deals 
with the formation of a relativistic bound state. We make the following ansatz: 
At some low value of — Qq, the transverse momenta of the partons may be 
ignored as a first approximation. The xb dependence of quark and anti-quark 
distributions at Qq, is then determined by solving a 1+1 dimensional model. In 
this model, quarks interact via a linear potential in the null coordinates. This 
is the simplest potential consistent with Lorentz covariance. 

We also suppose that the valence quarks carry a fraction / of the total 
baryon momentum at — Qq. The rest being carried by gluons, anti-quarks 
and sea-quarks. In a previous numerical study we determined the valence 
quark distribution in this interacting parton model. Here, we extend this model 
to include anti-quarks and also obtain analytic formulae for the valence quark 
distributions, within a variational approach. Based on the two parameters Qq 
and /, we predict the initial xb dependence of the iso-spin averaged quark and 
anti-quark distributions. Finally, to compare with experimental data at higher 
Q^, we use the solutions of this model as initial conditions for the DGLAP 
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evolution equations. 

We do not derive this model from a more basic theory. It is proposed merely 
as a phenomenological parton model for an approximate description of Deep 
Inelastic Scattering. 

Let us now give a brief introduction to our analysis of the interacting parton 
model. The partons are assumed to be relativistic particles interacting through 
a linear potential. The number of colors Nc is kept variable and we work mostly 
in the limit of a large number if colors. To simplify this many body problem 
we ignored the anti-quark degrees of freedom in |3( . The baryon wave function 
was determined by the principle that it minimizes the (mass)'^ of the baryon. 
Within a Hartree approximation, the valence quark wavefunction is the solution 
of a non-linear integral equation which was solved numerically in j3|. 

In this letter we first show that the true minimum of [mass)^ occurs for 
a configuration that consists only of valence quarks, in the chiral and large 
limits. The deviation of the anti-quark distribution from zero is measured by the 
dimensionless parameter ^ in addition to — corrections, which we estimate. 
Here m is the current quark mass and g a coupling constant. 

In order to determine the anti-quark content of the baryon, we perform a 
unitary transformation on the Fermionic Fock space, starting from a purely 
valence state. This transformation is like a Bogoliubov transformation which 
mixes positive and negative momentum states. It is sufficient to consider Bo- 
goliubov transformations parametrized by a single angle 6 which is determined 
by minimizing the [mass)'^ of the baryon. We find that in the large limit, 

2 

6 vanishes for zero current quark mass. For physically reasonable values of 2^, 
the anti-quarks carry less than a percent of baryon momentum. This is at an 
initial value of = Qq, at which transverse momenta are neglected. The 
isospin averaged quark and anti-quark distributions are found within a varia- 
tional approximation. 

To compare our results with experimental data, we evolve the distributions 
to higher values of via the DGLAP equations T. However, the gluon dis- 
tributions, which we have not determined, appear in the evolution equations. 
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It turns out that the difference between quarlc and anti-quark distributions 
(valence quark distribution: q^{xB,Q'^) = Hi3=u.Al'^ b ,Q^) - q^{xB,Q'^))) 
evolves independently of the gluon distribution to leading order. Moreover, if 
we ignore certain correlations, this difference is the average of the structure 
function measured in neutrino and anti-neutrino Deep Inelastic Scattering 
12 • Given the xb dependence of the parton distribution functions (PDFs) at 
an initial Qq, their evolution is determined by the DGLAP equation with 
the splitting function Pqg calculated perturbatively: 

dq^(xB,t) _ asjt) rl dy„V(„. j-\p ( '£b_\ 
dt ~ 2-nr JxB \y:'')^qq\ y )■ 

Here t = log{Q^/Q^). The normalization v{Q-^) = dxsq^ (xb^Qq) is deter- 
mined by integrating the DGLAP equation with initial condition v{co) = Nc = 3 
from = oo to Qq. Due to the large range involved, we determine i^iQa) 
to high order. Within our approximations, i'((5o) = Jq dxBF3{xB,Qa)i "^^^ch 
is given by the GLS sum rule 4 . If we denote the isospin averaged valence 
quark probabihty density as V{xb,Qo), then q^(x_B,Qo) = '^{Qo) V{xb,Qq)- 
In section 4 we compare our predictions for xF3{x, Q^) with experimental mea- 
surements by the CDHS and CCFR collaborations. Our predictions agree well 
with data for a choice of parameters I = \ and Qq — 0.4 GeV'^ for which 
i/((5q) = 2.25. This choice of parameters is consistent with phenomenological 
fits to data [21in|- However, it would be useful to know what the 'best-fit' values 
of these parameters are. 

An impressive discretized light-cone (DLCQ) analysis of 2 dimensional QCD 
was done by Hornbostel et. al. 6 . Our phenomenological parton model pro- 
vides a complementary physical approach to their more direct numerical diag- 
onalization of the hamiltonian. The 2-dimensional valence quark wave function 
we find reduces precisely to the one obtained in ^6. when we set f = 1. We do 
not find a similar concordance with the more conventional lattice QCD methods: 
the lightcone methods seem to incorporate the physical phenomena much more 
directly. Other DLCQ 7 calculations study the meson and glueball spectra of 
2d models. We focus on the baryon. For other approaches see for instance 
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2 Valence Parton Model 



Let us begin by reviewing the valence quark approximation. Rcf. 3 may be 
consulted for details. We assume that the momenta of the partons in the 
direction are large compared to the transverse momenta; which we ignore. We 
use null co-ordinates where the null momentum p — pq — pi is the basic 
kinematic variable. Then the kinetic energy of a free particle of mass m is 
Po = ^iP~^ ■^)- wave function of a quark will vanish for negative p while 

that of an anti-quark vanishes for positive p. 

If we ignore anti-quarks as in 3 , then the baryon wavefunction 
depends on the colors, flavours (a^ = 1, • • • , M) and null momenta of the Nc 
valence quarks. ip{{vi, ai,pi}) = e^^^...^^^^il;{{ai,pi}) since the baryon is a color 
singlet. Moreover, since the null momenta are positive, the sum of quark mo- 
menta cannot exceed the total baryon momentum P. In particular, the wave 
function must vanish for pi > P. Since the e tensor is anti-symmetric in color, 
the wavefunction must be symmetric in the remaining variables: partons behave 
like bosons in the momentum, spin and flavour variables. The ground state wave 
function is determined by minimizing the total energy 

Here /x^. — . — ^ is the effective mass of the parton 3 , which avoids a 
potential infrared divergence in the potential energy. Also, = . The 
simplest potential consistent with Lorentz invariance is linear, v(x) = which 
is also favoured by phenomenology g is a coupling constant with the 

dimensions of mass. Our predictions turn out to be independent of g in the 
chiral limit. In the ground state, we expect the Hartree ansatz 

= 27r(5(P - Y.^P^) Vu=\ i^i^i^Pi)- 

to be a good approximation. The valence quark wave function is normalized to 
have unit length, X^qLi Jq IV'Iq^jP)!^^^ — 1- It must also satisfy the momentum 

*See appendix to 9 for kinematics. 
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sum rule: lo PlV'(-P)l^|f = f^- f ^'^''^ ^^e fraction of baryon 

momentum carried by the valence quarks, which is roughly a half at low 
0. In these two formulae, we are ignoring correlations that are suppressed for 
large- A^'c. They differ from the exact formulae in the same way as the canonical 
ensemble differs from the micro-canonical ensemble in statistical mechanics. 

Since we are interested in the isospin averaged distributions, we will average 
over spin-flavour degrees of freedom. Thus we look for a wave function that is 
non-zero only for a single value of a: 'ip{a,p) — Sa^itp{p)- 

In second quantized language, this corresponds to the valence state \V >= 

• • • a^^^lO >. Here, a^^ creates a quark with color j in the state ^p. These op- 
erators satisfy canonical anti-commutation relations (CAR): {aiu,al'^} — 6f < 
u,v > with respect to the Dirac vacuum |0 > where all negative energy states 
are filled and positive ones empty: a^^ |0 >= and aj^^\0 >= 0. ip-{p) van- 
ishes for p > and i/j+ip) for p < 0. The Pauli principle requres that the 
density matrix pv{p,q) V \ -^a^\p)ai{q)\V > is a hermitean projection 
operator: J°^^pv{p,r) pv{r,q) |^ = pvip,q)- The eigenvalues of the den- 
sity matrix are the occupation numbers of particles, for a projection operator 
these are or 1 as required by the Pauli principle. For a state containing 
one baryon, the normal ordered trace of the density matrix is equal to one : 
tr{pY{p,q) + ^6{p,q){sgn{p) — 1)) = 1 by a use of the CAR. S{p,q) is the 
identity matrix. The above Hartree ansatz, — 'ip{p)ip*{q) + ^S{p,q){l~ 

sgn p) satisfies these constraints. 

A Lorentz invariant formulation is to minimize the mass of the Baryon 
per quark: 



-2 lo l-MP)\Hf + \ Poo dxdy\i:ix)my)\'^,\- - y\ 
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2.1 Analytic results in the large Nc limit 

In [3j we solved the integral equation for the minimization of energy numerically. 

There is in fact an analytic solution in the chiral > 0) and large Nc limits. 

The boundary condition is that ipip) must vanish for p > P. However, P is an 
extensive variable, P ~ Nc- So for Nc = oo, the valence quark wave function is 
not required to vanish for any finite value of p. If we use the intensive quantity 
P — , the analog of momentum fraction \s xb — , but the wave function is 
not required to vanish beyond = 1 . In order to compare directly with a wave 
function computed for Nc ~ 3, we pick the reference frame in which P — ^. 
The momentum sum rule then becomes /o°°p|'0(?')p5f — fP- It can be checked 
explicitly that V'(p) = \J jp^^^ is an exact solution to the integral equation 
for the minimization of baryon [mass)'^. Alternatively, we can calculate 
for this wavefunction and see that it is zero. The potential and self energies 
cancel each other. Thus, in this limit the minimum of (mass)^ actually occurs 
for a purely valence quark configuration. Since g ^ ^qcd ^ 200 MeV and 
Trim ™d ^5 — 8 MeV, this should be a good approximation provided the 
corrections are small. This is indeed the case, as we show below. For f — \, 
the valence quark density normalized to one is V{xb) — 6e~^^^. 



2.2 Leading order correction 

The leading order effect of finite Nc is to restrict the range of quark momenta 



to p < P. Now, (1 - ^) 



p_\n 



as n — *■ oo. Therefore, ip{p) = Cp"(l — 



p_\b 



< p < P should be a good ansatz for the ground state wave function for finite 
Nc- C is determined by normalization. The momentum sum rule implies that 
b — ~ 1 + ^("^ ^ The minimization of energy implies that a satisfies the 
transcendental equation 



7rm_ 



1 



+ Jo y- 



{i + yr + ii^yr 



oo dy 



which we derived in 131. In the limit of chiral symmetry, a — + 0. If valence 
quarks carry all the momentum of the baryon, / = 1, and V{xb) = [Nc — 
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Figure 1: Comparison of valence quark distribution xbQ^ {xb, Qq) obtained in 
the large- A'^c limit (thin curve) with the variational estimate xbQ^ {xb,Qo) after 
taking into account the leading correction (thick curve). 

1)(1 — xb)^"^'^, which is identical to the result obtained from DLCQ, reported 
in [S]. However, valence quarks carry only about half the baryon momentum, 
so that for Nc = 3 our variational estimate for the valence quark density is 
V{xb, Qo) — ^{^ ^ xb)*; this agrees well with our numerical solution from Ref. 

It is evident from Fig. 1 that the primary effect of the correction is to 
make the distribution vanish beyond -p = 1. The actual shape of the distribution 
is already well captured by our analytic solution in the large- A^c limit. 

3 Extension of Parton model to include Anti- 
quarks 

Now we turn to the anti-quark content of the baryon. As we will see below, the 
minimum of Baryon {mass)'^ occurs for a small anti-quark content when the 

2 

parameter Sj- is small. In order to determine it, we need to allow for states 
with negative momenta. However, only the difference between the quark and 
anti-quark numbers is conserved in the full theory. Therefore, the baryon must 
be in a linear superposition of states containing 77 anti-quarks and Nc + rj quarks, 
for 77 = 0, 1, • • • , 00. The energy of a state containing rj anti-quarks is 
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E., Ef=f'' I-p hiPi + ^Mr,i{'^j,Pj})\'Tikdpk 

where V'r;({^i>ft}) the wave function of such a state. {vi,Pi\^J:i^ refer to the 
colors and momenta of the quarks, while refer to anti-quarks. 

The wave function vanishes for negative quark momenta and positive anti-quark 
momenta. The total energy is the sum of the energies for each value of rj. 

Now we shall argue that the ground state of the baryon is to a good approx- 
imation determined by three orthonormal one-parton states ijj and which 
describe quarks and if)- which describes anti-quarks. We will continue to work 
in a factorized Hartree approximation, ignoring correlations except when they 
are required by the Pauli principle or color invariance. 

In the absence of anti-quarks, there are just Nc quarks whose wavefunction is 
completely anti-symmetric in color. Within the Hartree approximation, they all 
occupy the same single parton positive momentum state i/', which minimizes the 
{massY of the baryon. Now consider adding a quark and an anti-quark (97 = 1). 
The anti-quark will occupy the negative momentum state that minimizes 
the energy. Suppose all the Nc + I quarks occupy the state ip. Then, by the 
Pauli exclusion principle, the color part of the quark wave function must be 
totally anti-symmetric. However, there is no completely anti-symmetric tensor 
in A^c + 1 indices transforming under SU{Nc). Therefore, we are forced to 
introduce a new positive momentum state that must be orthogonal to the 
filled state ip. In order that this cost a minimal amount of energy, we expect 
tp+{p) to have only one more node than ip{p). As long as 77 < A^c) in the ground 
state, the quarks and anti-quarks occupy the states ■;/), i/i-. If 77 > Nc, we 
would have to introduce another pair of orthonormal states. However, we find 
that these additional corrections are very small. 

Now we shall argue that the configuration containing valence, sea and anti- 
quarks (say jT^S'A >) can be obtained by a unitary transformation acting on 
the valence quark state |y >: a Bogoliubov transformation. It must be unitary 
in order that < VSA\VSA >= 1. The operator that creates a quark in state 
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tjj+ and an antiquark in state V- is aj^_a^^^\ we sum over color indices to 
produce a color invariant state. Thus the unitary transformation we seek is the 
identity except in the two dimensional subspacc spanned by ■(/'+ and . Thus 
our variational ansatz is \VSA >= e^'"^'''-'**+ >. The density matrix of 

quarks in the new state can now be calculated: 

Pvsa{p, q) = '4'{p)i>{'l) - sin^ obP-{p)'4>-iq) - i>+{p)'>P+{Q)] 

-^sm29[tp^{p)tp+{q) + i:+{p)ip^{q)] + ^5{p,q){l - sgnp). 

Physical quantities are expressed most simply in terms of the normal ordered 
density matrix: M{p,q) = ^2p{p,q) + S{p,q){l — sgn p). For example, the 
baryon number density in momentum space is 

M{p,p)-M{-p,-p) = + sin2e - . 

This confirms the interpretation of ip- as the anti-quark wavefunction. The 
{massY is given by 



-\ /5p + 'i JZ, dxdy\M{x, yr\\x - y\ 

The variational quantities tp,'ip+,'tjj-, and 6 are determined by minimizing 
the baryon (mass)'^. In the ground state, we expect and 'tp-{p) to have 
no nodes (except possibly at the boundaries p = 0,P), while 'tp^ must have one 
more node. We estimate them variationally. 



3.1 Large Nc analysis 

Working in the Nc — > oo limit, the form of the analytic solution suggests the 



choice V'(p) = ^ ( f 



'i,V.+(p) = c+ f 



(!)• 


2 - Ci 




a ^ 



for p > Q and 



tp-ip) = '4'{~^p) p < 0. (For other ranges of p these functions must van- 
ish.) Ci is determined by the orthogonality condition while C, C+ are fixed by 
the normalization conditions. The variational parameter b determines the refer- 
ence frame. The Lorcntz invariant quantity A4^ is independent of b. Thus the 
variational principle will determine a and 9 and hence the wavefunctions. The 
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actual minimization of A^^ is a lengthy but straightforward calculation. Most 
of the energy integrals can be evaluated analytically and we do them using the 

2 

symbolic package Mathematica. We find that 0, a ^ as ^ ^ recovering 
the purely valence exponential solution. For a physically reasonable value of 
^ -^i^ 0.001, we estimate 9 = 0.02 and a = 0.035. This corresponds to 

S '^'qcd 

a small but non-vanishing anti-quark content in the baryon. 
3.2 Leading order correction 

As in the valence quark case, the leading effect is to make these wave func- 
tions vanish beyond p — P. We estimate this correction using the ansatz 
tjj{p) = Dp'^il - ]5)^'^/i+(p) = D+p'^ip - Di){l - pf for 1 > p > and 
ip-{p) — p) for — 1 < p < 0. Here P — 1 and for other ranges of p these 
functions must vanish. Di is determined by the orthogonality condition while 
D, arc fixed by the normalization conditions. For the choice Qq = OAGeV^, 
f = .001, we get = 0.02, a = 0.035, b = 2.175 for the variational 

parameters. The valence quark distribution is normalized to v{Qq) = 2.25 while 
the normalization of the anti-quark distribution is determined as a consequence 
to be i'{Qq) sin^ 6. Since sin^ 6 ^ 10^**, the primordial anti-quarks carry only 
about .01% of the baryon momentum. 

These results are identical to what we obtained from a more field theoretic 
point of view in [51^]. They also agree with the DLCQ analysis of [S] as pointed 
out in Section 2.2. However, the parton model point of view presented here is 
much simpler. Moreover, the GRV collaboration have obtained a reasonably 
good fit to Deep Inelastic Data for xb > 10^^ starting with a vanishing anti- 
quark distribution at an initial Qq ^ 0.2 GeV^. Our approximations are not 
expected to be valid for extremely low values of the momentum fraction, where 
the assumption that longitudinal momenta dominate, becomes questionable. 

Thus we find that the valence quark picture is quite accurate: the 'primor- 
dial' anti-quark distribution is very small. The anti-quark content is zero not 
only in the non-relativistic limit m » g but also (somewhat surprsingly) in 
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Figure 2: Comparison of predicted xFy, at = 13 GeV^ (solid curve) with 
measurements by CCFR (★) at 12.6 GeV^ and CDHS (0) at 12.05 < < 14.3 
GeV2. Ql = 0.4 GeV2 and / = i. 

the chiral hmit to = when Nc oo, with corrections being smalL Never- 
theless, a substantial anti-quark content is generated by evolution. 

4 Comparison with experimental data 

Finally, we compare with experimental data at higher values of . The DGLAP 
evolution equation is integrated numerically. We set A^c = 3, Aqcd = 200 
MeV and the current quark mass to = 0. The parameters /, Qq should be deter- 
mined by a best fit to experimental data. For now we assign to them reasonable 
values / = 5 and Qq — 0.4 GeV^ which are consistent with phenomenological 
fits to data. From the GLS sum rule we get j/(0.4 GeV^) = 2.25 g]. In Fig. 2, 
we show a comparison with xF^ measurements by the CDHS and CCFR col- 
laborations [12 at Q2 _ 13 GeV^. The small range 0.4 < < 13 GeV^ over 
which we are evolving justifies the use of the leading order DGLAP equation. 
The plot shows that our prediction agrees reasonably well with the experimental 
measurements. 
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